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Abstract. - Our interest goes to the power injected in a heated granular gas and to the 
possibility to interpret it in terms of entropy flow. We numerically determine the distribution 
of the injected power by means of Monte-Carlo simulations. Then, we provide a kinetic theory 
approach to the computation of such a distribution function. Finally, after showing why the 
injected power does not satisfy a Fluctuation Relation a la Gallavotti-Cohen, we put forward 
a new quantity which does fulfill such a relation, and is not only applicable in a variety of 
frameworks outside the granular world, but also experimentally accessible. 



In the realm of non-equilibrium statistical physics, universal or generic results are scarce. 
Landmarks in the field are Einstein's seminal contribution [1] exhibiting a relation between 
particle current fluctuations and a response function, followed by Onsager's reciprocity [2] 
and the Green-Kubo relations [3,4]. The recent discovery by Evans, Cohen and Morris [5] 
formalized into a mathematical theorem of general scope by Gallavotti and Cohen [6] of a 
symmetry property, that we will call Gallavotti-Cohen relation (GC), bearing on the entropy 
flux distribution therefore stands as a major progress. Similar results have been established 
for Markov processes by Kurchan [7] and Lebowitz and Spohn [8], with much lighter math- 
ematics than for dynamical systems. However, mathematical as well as physical difficulties, 
that we shall later describe, have prevented this result to find its way towards experimental 
confirmation. 

From the point of view of probing the validity of the Gallavotti-Cohen theorem outside 
its mathematical domain of validity, granular gases will prove instrumental. This was indeed 
already realized by Aumaitre et al. [9]. Yet, in the absence of a first-principle definition of 
entropy, a physical one, albeit heuristic, was proposed : the energy injected by the thermostat 
divided by the granular temperature. This is precisely the quantity simulated in [9], and 
considered in a recent experimental work [10]. Besides, the idea of using a macroscopic 
and global observable for characterizing the state of a system, instead of resorting to local 
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probes (velocity field, correlations, structure factor) has proved a valuable tool per se for the 
identification of generic features in non-equilibrium systems [11]. 

In this letter we consider a paradigmatic model of a granular gas maintained in a steady- 
state by external heating, namely hard-sphere (or hard-disk) particles undergoing inelastic 
collisions, each of these particles being independently subjected to a random force Fj which 
we take to be Gaussian distributed with variance T. The equation of motion for particle i 
with velocity Vj then reads dvi/dt = F co n + Fj where F co u is the force due to inter-particle 
collisions, and acts at contact only [12]. This model preserves essential features of standard 
experimental setups, like the inelastic collisions, and a heating mechanism independent of the 
particles' velocities. It moreover bypasses experimental difficulties like the lack of transla- 
tional invariance and has been the subject of intense investigations both on the numerical and 
analytical sides [12], therefore its putative limitations are controlled. 

The results we have obtained arc as follows: First, by means of Monte-Carlo simulations, 
we pinpoint the numerical hazards that pave the way to the full determination of the injected 
power probability distribution function (pdf) and its symmetry properties. Second, we show 
how kinetic theory can be extended to describe the large deviation function of the injected 
power thus allowing for an explicit test of GC. This constitutes the first analytical result 
of a distribution function for a global observable in a many-body, non Gaussian, far from 
equilibrium system. Finally, we explain the reasons why the injected power cannot satisfy a 
Gallavotti-Cohen relation. This leads us to proposing a new quantity, with the properties of 
a Gibbs entropy flow, that not only avoids the aforementioned difficulties, but that should 
also be accessible in specific granular gases experiments, and could possibly be generalized to 
other non-equilibrium systems. 

We define the time-integrated injected power over the time interval [0, t] as the total work 
W(t) provided by the thermostat: 

W(t)= [ dt VF 4 -v 4 . (1) 
Jo 

The granular temperature is defined by T g — 1//3 — (wf)/d where d is the space dimension 
and the angular brackets denote a statistical average in the steady-state. The central object 
of our study is the pdf of W(t), P(W, t), and its related large deviation function 

TToo(u') = lim ir t (w) with n t (w) — — In P{wt 1 1). (2) 

t— >oo t 

The generating function P(A,t) = (e~ AW ) and its related large deviation function yu(A) = 
lim^oo j In P(A, t) will provide equivalent information, as we may usually go from one to the 
other by Legendre transform, ^^(w) = max,\{^i(A) + Xw}. 

We now recall the content of the GC relation. In phcnomcnological thermodynamics [13], 
one writes an evolution equation for the entropy 5, 

^ = t- [dVV-Js (3) 

^ J system 

which, in a steady-state, expresses a balance of an entropy production term (a measure of 
intrinsic irreversibility) with an entropy flux 3$ produced by an external drive (typically, a 
boundary condition imposing energy or matter flow). The central result of Gallavotti and 
Cohen is to have shown in a particular setting that the pdf P(S, t) of a quantity S(t) playing 
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the role of the time- integrated flux term J dtdT^V ■ Js verifies, in the limit of asymptotically 
large times, 

P(-S,t) 1 ' 

Nevertheless, the assumptions underlying GC are microscopic reversibility and the identifica- 
tion of S with the phase space contraction rate. In granular materials, the former is violated 
and the latter definition is doomed to fail, since the phase space volume can only decrease 
if the heating mechanism is velocity independent. This has led to the idea that /3W could 
possibly play the role of S [9] . In the heated granular gas the fluctuating total kinetic energy 
E{t) = J2i v i 2 /2 varies according to 

E(t) - E{0) = W(i) - V{t) (5) 

where T>(t) > is the energy dissipated through collisions up until time t. Of course on 
average (W(t)) = (T>(t)) — 2dTt > 0, yet there will be individual realizations for which W 
will be negative. If the GC relation Q held for W, it would take the following form [14] 

7Too(w) - 7roo(-w) = f3w, ^i(A) = n([3 - A) (6) 

Remark, though, that the left-hand side in |5j) is bounded in time (while W and T> grow 
typically linearly with time), therefore [15] both the injected power and the dissipated power 
have the same (cumulant) generating function. A straightforward consequence of that fact 
is the absence of w < events at large times: ^^(w < 0) = — oo, hence © cannot be 
correct [16]. 

Let us see now how this rigorous fact comes about in numerical studies. We have simulated 
a model of N inelastic hard disks under the effect of a homogeneous thermostat. The collisions 
between two disks i and j conserve the total momentum and reduce the normal component 
of the relative velocity i.e. (v^ — v'-) • n = — a(vj — Vj) • n where the primes mark the post- 
collisional velocities and n is the direction joining the colliding particles. The driving force 
F t mentioned above obeys {F?(t)F$(tf)) = 2T5 lj S' lS S(t - t'). Molecular Chaos is enforced 
through the use of a Direct Simulation Monte Carlo [17] algorithm. The restitution coefficient 
a takes values between and 1 (elastic gas). The gas rapidly reaches a stationary state 
with a granular temperature T g — AT /(I — a 2 ), where the mean free time is used as a time 
unit, which will be the case in the subsequent analysis. We have measured the fluctuations 
of the work exerted by the thermostat W(t) as defined in Eq. using an integration time 
t. P(W, t) clearly shows non-Gaussian tails, which can be regarded as a hint that the true 
large deviations are being probed (since central limit theorem requires small deviations to be 
Gaussian). The resulting function TVt(w) is shown (see Fig. to verify almost perfectly a 
relation such as When decreasing a a slight departure (about 10%) from the slope /3 can 
be recognized. And yet the plot of TT t (w) — ir t (—w) stays linear in w. The inset of Fig. ^ is of 
particular interest: it is numerically the proof that the slope does not depend upon the time of 
integration for the considered range of t values. More interestingly, it highlights the dramatic 
decay of the number of observable negative events as t is increased: even with such a simple 
and fast to simulate model, obtaining large statistics at large t is difficult. This observation 
poses a crucial question: has the numerical investigation reached the true asymptotics, i.e. 
can we believe that n t (w) ~ ^^(w) for the values of t shown in the figure? In the following 
we exploit further numerical studies and a novel analytical approach to demonstrate that a 
time much longer should be waited and that this first numerical result is misleading. It is 
nevertheless important to discover that at values of t small enough to observe some negative 
injected power the GC relation appears to be satisfied. 
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Fig. 1 - Color online: Plot of itt(w) — n t (—w) as a function of /3w. The dashed line has slope 1, 
the dotted line has slope 1.1. The inset contains the same graph for different values of t, for the case 
a = 0.9, N = 100 and T = 0.5. 

Fig. 2 - The large deviation function -k^, of the rescaled quantity W/ { W) obtained from the Legendre 
transform of the largest eigenvalue p(X) (see text and [18] for details). The inset shows the dimension- 
less third cumulant (VV 3 ) C = (W 3 ) c /3 2 /(VV) relaxing towards its predicted theoretical value, which 
is 8. The part of the curve close to w/(w) is magnified to show the agreement between simulations 
and theory on the limited accessible range of values at t — 40. 



We now sketch the analytical strategy that allowed us to find the pdf of the injected power. 
We start from an extended Liouville equation for the probability p(Tj^, W, t) that the system 
is in state Tjy with W(f ) = Wat time t. The second step is to convert the Liouville equation in 
terms of the generating function p(Tn, A, t) — J dWe _AW p. Recall that, for A = 0, projection 
of the Liouville equation onto the one-particle subspace yields an equation coupling the one- 
particle distribution function to the two-particle distribution function, which is factorized 
through the molecular-chaos hypothesis. In our case, the situation is very similar when it 
comes to determining the largest eigenvalue p(X) of the evolution operator for p (which, in 
physical words, is the generating function for the cumulants of W(t)/t); projecting onto the 
one-particle subspace yields an equation coupled to the two-particle subspace. By means of a 
molecular-chaos like closure procedure we arrive at an equation for both the eigenvalue p and 
its related eigenfunction. The extensive mathematical analysis of this program is cumbersome 
and will be reported elsewhere [18]. We concentrate here on the results that are as follows. 
First we find that the large deviation function of W = wt has the graph depicted in Fig. |21 
The tails are given by 

7r oo (w->0 + ) w" 1/3 , 7TCX, (to ->• +oo) w (7) 

with, as expected, no w < contribution. While the w — > + regime appears to be thermostat- 
dependent, the exponential right tail of P(W,t) seems to be a robust property related to the 
presence of a branch cut in the complex A plane for p(X). 

We now come back to our previous numerical results. The observation of a clean linear 
behaviour for 7Tt(iw) — 7Tt(— w), with slope (3, apparently independently of time (see inset of 
Fig. nj, as though GC were satisfied, is in contrast with our analytical results indicating 
that the truly asymptotic n^w) has no w < contribution. At the level of simulations or 
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experiments it is impossible to distinguish between non-Gaussian tails due to spurious short 
time effects and real asymptotic large deviation tails. Nevertheless the numerical study of 
cumulants with time, which are integrals and therefore benefit from a larger statistics, is 
decisive both to prove the validity of our analytical approach and to get rid of all doubts 
about which is the true asymptotic regime. We see (inset of fig. EJ that the third cumulant 
reaches a stable value at times of order ~ 50, much larger than the times used in Fig. ^ 
Moreover this stable value is in very good agreement with the value expected from our theory 
(W 3 ) c = — t [d 3 /.i(\)/d\ 3 ]\ = Q. We recall that this estimate of the third cumulant is highly 
non-trivial. At such large t values, however, the numerically accessible range of w is very 
small, limiting the possibility of a detailed comparison between analytics and numerics (see 



Having reached the conclusion that the injected power cannot fulfill the GC relation, in 
this final paragraph we propose an alternative quantity which is, by its very definition, a 
Gibbs entropy flow, and satisfies the GC fluctuation relation. The sequel applies both to 
the original experimental system in which particles contained in a closed box are vigorously 
shaken and to our random thermostat. We now tag one of these particles that we follow 
along its path, which is similar in spirit to the work reported in [19]. The rest of the particles 
act as a thermostat for the tagged particle. Note however that the states of the particle 
bath evolve according to dynamical rules that do not fulfill the detailed balance condition. 
Similarly, the trajectory in the phase space of the tagged particle does not fulfill the detailed 
balance condition (this is a key difference with [19]). It is now possible to view the time 
evolution of the tagged particle as a Markov process. In the tagged particle phase space, to 
every non-zero transition rate between two velocity states one can associate a non-zero rate 
for the backward move. This is the core of the difference with following the dynamics in the 
phase space of the whole system. There, due to irreversible microscopic dynamics (inelastic 
collisions), for a given forward- in-time trajectory in phase space, there is no corresponding, 
however extremely unlikely, backward-in-time trajectory. Hence for the tagged particle it is 
perfectly legitimate to define an integrated entropy flow a la Lebowitz and Spohn [8] (see also 
Maes [20] or Gaspard [21] for alternative presentations): 



where n(t) is the number of collisions undergone by the tagged particle over the time interval 
[0,i], Vi is its velocity after the i-th collision, and K(v — > v') is the velocity transition rate 
of a particle undergoing a collision. Expression (JHJ may be experimentally accessible in the 
following way: the monitoring of a large number of collisions (and of the velocities before 
and after collision) first gives an estimation of the velocity transition rates K(v — > v'), that 
can be tabulated with an appropriate discretization of velocities; it is then possible to follow 
one given tagged particle, to measure its successive pre- and post-collision velocities and thus 
to compute S. In the elastic limit for the thermostat the transition rates K(v — > v') verify 
the detailed balance condition with respect to a Maxwellian tagged particle velocity pdf [19], 

hence there is no entropy flow anymore, and S(t)/t (the tagged particle is in equilibrium 
within its own phase space). Fully explicit expressions for the transition rate K(v — > v'), and 
hence for S(t) [18] can be derived for an inelastic thermostat, e.g. using Sonine expansions. 
By construction, the pdf of the instantaneous entropy flow ijSJ will verify a GC relation, with 
the ensuing consequences in terms of Green-Kubo relations. 

We have provided the first computation of the large deviation function of the injected 
power in a system driven far from equilibrium, for which no general theory has hitherto 
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been coined, as opposed to systems weakly driven out of equilibrium [22]. This computation 
allowed us to pinpoint numerical limitations and yet to extract reliable data. We have shown 
why the injected power cannot satisfy a GC relation, thus revealing a scenario sharing many 
features in common with the equilibrium toy model considered by Farago [23]. However, even 
if 7TOO, the large t limit of nt(w), cannot exhibit a GC symmetry due to the absence of a 
negative tail, it appears that from a practical point of view, %t(w) seems to obey a GC-like 
relation (see Fig^]). An important consequence of the analytical work outlined here is that 
such a "finite-time" property cannot be considered as an extension of GC theorem. Finally 
we have put forward a new approach, based upon a Lagrangian point of view, that leads 
to a definition of an entropy flow possessing all the properties requested by the systematic 
approach of Lebowitz and Spohn [8], and which may be experimentally accessible. We are 
confident that not only the new extensions of kinetic theory developed here, but also the 
proposal for an entropy function, will trigger new theoretical investigations and will lend 
themselves to experimental confirmations far beyond the granular community, in the spirit of 
recent experiments on turbulent flows resting on Lagrangian viewpoint investigations [24] . 
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